Response traits and stability in the context of a pulse perturbation. Are fundamental niche- or realised niche-based response traits more informative?

1 Introduction

In the context of a pulse perturbation we ask:

  • will fundamental niche-based response traits, such as those measured according to Ross et al (2023), explain variation in community stability?
  • is a realised niche-based response trait able to explain more variation in community stability?

We explore these questions using a simulation model of a multi-species model, in which species intrinsic growth rate is a bell-shaped function of temperature, and each species can have a different temperature optimum. A pulse perturbation is created by briefly setting temperature lower than the control temperature, and then returning it to control temperature.

We then manipulate the strength of inter-specific interactions and examine how that affects the explanatory power of the two types of response trait.

Predictions:

  • Realised and fundamental niche-based response traits will provide equal and high explanatory power when there is no or very weak interspecific interactions.
  • Increase in strength of interspecific interactions will decrease the explanatory power provided by fundamental niche-based response traits.
  • Increase in strength of interspecific interactions will not decrease, or will decrease relatively slowly (compared to fundamental niche-based response traits) the explanatory power provided by realised niche-based response traits.
  • Realised niche-based response traits will provide greater explanatory (than fundamental niche-based response traits) when interspecific interactions are relatively strong.

2 Methods

2.1 The multispecies model

We have \(S\) species that can be interacting and whose vital rates are temperature dependent. We assume density-dependent birth rate, (\(B\)), and death rate (\(D\)) in a discrete-time version of the classical Lotka–Volterra model [@de2013predicting; @vasseur2020impact] to get instantaneous growth rate, \(\tilde{r}_{i}(t)\), for species, \(i\), in year \(t\):

\[\begin{equation}\label{eq.r} \tilde{r}_{i}(t) = ln N_{i}(t+1) - ln N_{i}(t) = B(N_{i}(t),N_{j}(t),T(t)) - D(N_{i}(t),N_{j}(t),T(t))(\#eq:r) \end{equation}\]

Here, \(N(t)\) represents the biomass at year \(t\), and \(i,j\) are indices for two different species.

The per-capita birth and death rates for \(i^{th}\) species are represented as:

\[\begin{equation}\label{eq.B} B_{i} = b_{0,i}(T)-\beta (N_{i}+\sum_{i \neq j = 1}^{S} \alpha_{ij}N_{j})(\#eq:B) \end{equation}\] \[\begin{equation}\label{eq.D} D_{i} = d_{0,i}(T)-\delta (N_{i}+\sum_{i \neq j = 1}^{S} \alpha_{ij}N_{j})(\#eq:D) \end{equation}\]

where, \(\beta\) and \(\delta\) are density-dependent constants, \(\alpha_{i,j}\) is the competition coefficient between species \(i\) and \(j\), and

\[\begin{equation}\label{eq.b0} b_{0,i}(T) = a_{b} e^{-(T-b_{opt,i})^2/s_{i}} (\#eq:b0) \end{equation}\] \[\begin{equation}\label{eq.d0} d_{0,i}(T) = a_{d} e^{z_{i}T} (\#eq:d0) \end{equation}\]

with \(a_{b}\), \(a_{d}\) as intercepts, and for \(i^{th}\) species, \(b_{opt,i}\) is the temperature that optimizes birth rate, \(s_{i}\) governs the breadth of the birth function, \(z_{i}\) scales the effect of temperature (in °C) to mimic the Arrhenius relationship.

Substituting Eqs. (\(\ref{eq.B}\)) - (\(\ref{eq.d0}\)) in Eq. (\(\ref{eq.r}\)), we get the following,

\[\begin{equation}\label{eq.r2} \tilde{r}_{i}(t) = r_{m,i} \left( 1-\dfrac{\sum_{i, j = 1}^{S} \alpha_{ij}N_{j}}{K_{i}} \right) (\#eq:r2) \end{equation}\]

where, \(\alpha_{ii} = 1\), \(r_{m,i} = (b_{0,i} - d_{0,i})\) is the intrinsic (maximum) rate of increase and \(K_{i} = r_{m,i}/(\beta+\delta)\) is the carrying capacity for the \(i^{th}\) species, respectively.

If abundance drops very low then it is reset to 1. May need to change this, or at least present a good ecological justification.

2.2 Model parameters

The following need to be specified for a community (in addition to the number of species):

Property Level Notes
a_b Species Intercept of intrinsic growth rate - temperature function of species i.
b_{opt,i} Species Temperature at which intrinsic growth rate is highest for species i.
s_i Species Width of intrinsic growth rate - temperature function of species i.
a_d Species Intercept of death rate - temperature function of species i.
z_i Species Slope of death rate - temperature function of species i.
alpha_{i,j} Interaction Strength of interspecific effect of species j on species i.

2.3 Experimental treatments

2.3.1 Pulse pertubation treatment

The context of our entire study is a pulse perturbation. Temperature is constant in one treatment, and is pulse perturbed in another. For simplicity, we always make the perturbed temperature lower than the control temperature. We do not manipulate or change anything about the pulse perturbation. The graph below shows the control temperature in black and the pulse in red. Simulations are run for 10’000 time steps before the pulse is applied and run for another 1’000 after the pulse starts.

2.3.2 Creating a community

In order to create a community we define a set of \(S\) species and in the following only two things differ among the species:

  1. The optimum temperature for growth, such that some species have maximum intrinsic growth rate at higher temperatures and some have maximum intrinsic growth at lower temperatures. The optimum temperature for growth for each species in a community is drawn from a uniform distribution with a specific mean and range.
  2. The interspecific interactions in the interaction matrix \(alpha_{i,j}\). These are drawn from a normal distribution with mean zero and standard deviation \(sd(alpha_{i,j})\).

2.3.3 Community composition treatment

We create communities that vary in their composition by changing the mean of the uniform distribution from which species’ optimum temperature for growth are drawn. If we set the mean low, close to the temperature of the perturbation, then the perturbation as a positive effect on many of the species intrinsic growth rates. In contrast, if we set the mean close to the control temperature then the perturbation has mostly negative effects on the species’ intrinsic growth rates. We show some communities in the Results section of this report.

2.3.4 Strength of species interactions treatment

When we set \(sd(alpha_{i,j})\) = 0 then there are no interspecific interactions. Larger values of \(sd(alpha_{i,j})\) mean that interspecific interactions are stronger.

2.4 Analyses

2.4.1 Species response traits

We calculated two types of species response trait. First was a response trait based on the temperature response curve of a species (i.e., as in Ross et al (2023)). For this trait, we calculated the effect of the perturbation on intrinsic growth rate (perturbation temperature intrinsic growth rate - control temperature intrinsic growth rate). This is analogous to the calculating the first derivative (slope) of the temperature response curve.

Second was a response trait based on the observed response of the species to the pulse perturbation in the community context. Hence this response trait depends on the direct effect of the pulse on a species, and also the indirect effects via other species in the community.

2.4.2 Community response diversity

Currently we calculate:

  • IGR-effect trait mean, variance, dissimilarity response diversity, and divergence response diversity.
  • Species AUC trait mean and variance.

2.4.3 Community stability

This was calculated as the effect of the pulse perturbation on the total community biomass. This effect varied through time (e.g., often decreased monotonically after the pulse ended). The effect was calculated as the control - perturbed abundance and could conceivably be sometimes positive and sometimes negative. Examples of the dynamics of total biomass in control and perturbed treatments are given in the Appendix. A schematic figure is just below.

In one measure of community stability, we took the sum (over time) of the effect. If the effect was always positive, the summed effect would be positive and if the effect was always negative then the summed effect would be negative. A small (close to zero) summed effect could result from a balance of positive and negative effects, or just very small but consistent positive (or negative) effects. Stability is greatest when the summed effect is zero. More negative and more positive values imply lower stability. We call this raw community stability.

Another measure was the sum of the absolute effect size (OEV, Overall Ecological Vulnerability Pablo et al 2021. This will only be small (closer to zero) when the perturbation has a consistently small effect (positive, or negative, or mix of positive and negative). A mix of large positive and large positive effect would lead to a large value of this measure.

community stability schematic

3 Results

3.1 Explanatory power

$spp_RR_calc_threshold
[1] 0.01

$interactions
[1] "only_comp"
[1] "pack3"

And here we show how much variability in community stability is explained by the community mean of the species response traits, for realised niche-based response traits, and for fundamental niche-based response traits. How much variability explained is that of a GAM.

Seems surprising that the fundamental niche-based response traits always (on average) provide greater explanatory power. This contradicts expectations in the Introduction section (Section 1).

3.2 Response trait and community stability

In the following various graphs is the relationship between community mean of species response trait (realised or fundamental) (x-axis) and community stability (y-axis).

Each panel is for a different strength of interspecific interaction.

3.3 Response balance

Ross et al (2023) proposed that response diversity not just be measured by the spread of response trait values, but also by including if there are some negative and some positive response trait values. They made a measure, called divergence, that is zero when all response trait values are negative or all are positive, and is 1 when the most negative value is equal in magnitude to the most positive value.

In the results above we see that community stability is greatest when the mean of the perturbation effect on species intrinsic growth is zero. A mean of zero results from balance in the negative and positive response trait values, i.e., the sum of the positive values equals the sum of the negative values.

Hence, at least in the context of a pulse perturbation we see that the response balance, which is calculated as the mean of effects of the pulse on intrinsic growth rate, is an excellent predictor of community stability, at least when interspecific interactions are weak.

We can go a little further, and predict that response balance (i.e., mean_species_igr_effect = 0) will be greatest when the expected temperature optima of the species is half-way between the control and perturbed temperature. This is exactly what occurs in the data, as shown in the graph just below.

One set of simulations includes communities that have different values of the expected range of b_opt and show that this does not affect community stability to the pulse perturbation.

3.4 Some individual communities

Examples 1, 2, and 3 are with no interspecific interactions. The first and second are with relatively low variation in species responses, whereare the third has higher variation (half positive and half negative). Notice how the resistance of total biomass of the third community is greater than the resistance of the first two.

Examples 4 has relatively weak interspecific interaction, and example 5 has relatively strong interspecific interaction. The black boxes in some of the individual species responses are the addition of 1 abundance unit when the abundance drop too low (see last sentence in section Section 2.1).

3.4.1 Example 1

No inter-specific interactions (sd(alpha_ij) = 0) and low average optimum temperature (mean(b_opt) = 17), so the pulse has positive effect on intrinsic growth rate of most species.

species_id species_RR_AUC igr_pert_effect
Spp1 145.37993 0.2977747
Spp2 -24.83772 -0.1068233
Spp3 151.51334 0.2949134
Spp4 25.01589 0.0828334
Spp5 150.76799 0.2957067
Spp6 37.51282 0.1222027
Spp7 53.59615 0.1690748
Spp8 49.88283 0.1588482
Spp9 152.14160 0.2939789
Spp10 15.81453 0.0531594

3.4.2 Example 2

No inter-specific interactions (sd(alpha_ij) = 0) and high average optimum temperature (mean(b_opt) = 20), so the pulse has negative effect on intrinsic growth rate of most species.

species_id species_RR_AUC igr_pert_effect
Spp1 32.9589393 0.1080246
Spp2 -42.2535099 -0.2912424
Spp3 -39.4152687 -0.2380594
Spp4 -41.4311479 -0.2742790
Spp5 0.2774118 0.0009767
Spp6 -38.8734471 -0.2298022
Spp7 -41.3172661 -0.2719812
Spp8 19.9556752 0.0665788
Spp9 17.2907469 0.0579568
Spp10 -36.7273678 -0.2015994

3.4.3 Example 3

No inter-specific interactions (sd(alpha_ij) = 0) and intermediate average optimum temperature (mean(b_opt) = 19), so the pulse has negative effect on intrinsic growth rate of some species, and positive effect on others.

species_id species_RR_AUC igr_pert_effect
Spp1 76.83214 0.2224417
Spp2 -34.43342 -0.1771854
Spp3 -40.07936 -0.2489685
Spp4 -40.90983 -0.2640022
Spp5 17.61857 0.0590200
Spp6 -40.94920 -0.2647559
Spp7 57.04435 0.1781744
Spp8 -41.80079 -0.2818977
Spp9 -17.16628 -0.0676524
Spp10 -40.47937 -0.2560059

3.4.4 Example 4

“Weak” inter-specific interactions (sd(alpha_ij) = 0.1) and intermediate average optimum temperature (mean(b_opt) = 19), so the pulse has negative effect on intrinsic growth rate of some species, and positive effect on others.

species_id species_RR_AUC igr_pert_effect
Spp1 NA 0.1701316
Spp2 -13.767747 -0.0708877
Spp3 188.781686 0.1147783
Spp4 -10.057010 -0.0539152
Spp5 129.983516 0.1797482
Spp6 -44.880329 -0.2328473
Spp7 -28.179341 -0.1048868
Spp8 -18.485129 -0.1147155
Spp9 6.526969 -0.0330307
Spp10 -41.989491 -0.2003556

3.4.5 Example 5

“Strong* inter-specific interactions (sd(alpha_ij) = 0.3) and intermediate average optimum temperature (mean(b_opt) = 19), so the pulse has negative effect on intrinsic growth rate of some species, and positive effect on others.

species_id species_RR_AUC igr_pert_effect
Spp1 -35.48525 -0.1687307
Spp2 -15.65081 -0.0777735
Spp3 NA 0.0066435
Spp4 NA 0.1961017
Spp5 NA 0.2014327
Spp6 -42.24509 -0.2618709
Spp7 NA 0.2208124
Spp8 NA 0.1766150
Spp9 -23.84014 -0.1023592
Spp10 NA 0.1786434

3.4.6 Code for checking a specific case

Comm-1241-rep-1

species_id species_RR_AUC igr_pert_effect
Spp1 1.488498 -0.0278086
Spp2 NA -0.0373169
Spp3 NA -0.0352354
Spp4 5.975944 -0.0696115
Spp5 -36.936201 -0.1619554
Spp6 5.330759 -0.0602790
Spp7 -48.856549 -0.2806594
Spp8 -40.313351 -0.2411587
Spp9 -20.230398 -0.1185603
Spp10 22.328828 -0.0924903

4 Appendix

4.1 Check various relationships

5 From Charly